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Abstract: 

Extensive studies on thin films indicated a generic cubic current-voltage I — V dependence as a salient feature of charge transport 
by tunneling. A quick glance at / — V data for molecular junctions suggests a qualitatively similar behavior. This would render 
model-based studies almost irrelevant, since, whatever the model, its parameters can always be adjusted to ht symmetric (asym¬ 
metric) I — V curves characterized by two (three) expansion coefficients. Here, we systematically examine popular models based 
on tunneling barrier or tight-binding pictures and demonstrate that, for a quantitative description at biases of interest (V slightly 
higher than the transition voltage V)), cubic expansions do not suffice. A detailed collection of analytical formulae as well as their 
conditions of applicability are presented to facilitate experimentalists colleagues to process and interpret their experimental data 
by obtained by measuring currents in molecular junctions. We discuss in detail the limits of applicability of the various models 
and emphasize that uncritically adjusting model parameters to experiment may be unjustihed because the values deduced in this 
way may fall in ranges rendering a specihc model invalid or incompatible to ab initio estimates. We exemplify with the bench¬ 
mark case of oligophenylene-basedjunctions, for which results of ab initio quantum chemical calculations are also reported. As 
a specihc issue, we address the impact of the spatial potential prohle and show that it is not notable up to biases V ^V), unlike at 
higher biases, where it may be responsible for negative differential resistance effects. 


Keywords: molecular electronics; electron transport; single-molecule junctions; tunneling barrier; tight binding models; transi¬ 
tion voltage spectroscopy 


1 Introduction and background 


The roughly parabolic shape of the conductance GiV) = 
dl/dV, or the related cubic dependence of the current (I) on 
bias (V), was considered a prominent characteristic of trans¬ 
port via tunneling. This conclusion emerged from extensive 
studies on a variety of macroscopic thin him junctions of 
oxides, insulators, superconductors up to relatively large bi- 
ases.itt^ A quick glance at I — V measurements in a variety 
of molecular junctions may convey the impression that this 
cubic dependence is satisfactory also for such systems.— As il¬ 
lustration, we have chosen in Fig.[TJi a measured I — V curve,^iS 
for which htting with a cubic polynomial looks particularly accurate. 
The view based of such third-order Taylor expansions (cf. eqn O) 
was able to qualitatively describe a series of interesting aspects re- 
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lated to charge transport in molecular junctionsi^ 

liv) = h{v) + ff{y^y,h{v) = G{y + c2V^^ 

GiV) = ^/(V') = G(l+3c2l/2)+^(v'4j (1) 

Above, G = Ivoiv^qI/V is the low bias conductance.— 

However, if the approximation of eqn O were quantitatively 
adequate, it would be deceptive for a model-based description of 
transport in molecular junctions. Resorting to simple phenomeno¬ 
logical models enabling expedient data processing and interpreta¬ 
tion of charge transport through nanoscale devices represents com¬ 
mon practice in molecular electronics. The validity of a certain 
model/mechanism for transport in a given system/device is often as¬ 
sessed by considering its ability to ht the measured current-voltage 
I — V curves. Adopting this “pragmatic” standpoint, the problem en¬ 
countered is that an approximate generic parabolic shape of the con¬ 
ductance could be inferred at not too high biases regardless the tun¬ 
neling model.^ Whether the electron wave function tunnels across 
a structureless average medium modeled as a tunneling barrier or 
through tails of densities of states of one (frontier) or a few off- 
resonant molecular orbital levels, one could “appropriately” adjust 
a few parameters, and the third-order expansion of the current for¬ 
mula I = I{V), often available in closed analytical forms from lit- 
erature, 4, 6 , 11,12 could relatively easy be “made” to ht the usually 
featureless I — V curves measured. Then the quality of the htting 
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alone cannot be invoked in favor of one tunneling mechanism out of 
other possible mechanisms underlying the various phenomenological 
models. 

The inspection of Fig. & may indeed convey the impression of an 
overall very good agreement between experiment and eqn Q- How¬ 
ever, a more careful analysis reveals that the range of the highest 
biases that can be sampled experimentally is quantitatively not so 
satisfactorily described. Typically, these V-values are only slightly 
larger than the transition voltage Vf.— The difference in V/-values 
extracted from the experimental and fitting curves of Fig.lJJ) exceeds 
typical V,-experimental errors. 

The range V ~ V, is of special interest. As discussed in Section 
[ 2 I Vt represents an intrinsic property of a junction out of equilibrium; 
it characterizes the (high) bias range exhibiting significant nonlinear¬ 
ity. Therefore, even more than an overall fitting of / — V transport 
data, it is important for theory to correctly understand/reproduce V). 
And, as illustrated in Fig. a and c, to this aim, the theoretical de¬ 
scription should go beyond the cubic expansion framework provided 
by eqn Q. Emphasizing this aspect in the analysis of the various 
transport models considered below, which are among the most popu¬ 
lar in the molecular electronics community,— represents 
one main aim of the present study. 

Emerging from idealized description of reality, the models utilized 
in transport studies, like those used elsewhere, inherently have a lim¬ 
ited validity. The various model parameter values are subject to spe¬ 
cific restrictions, which represent intrinsic limitations for the model 
in question; merely succeeding to provide good quality data fitting 
is meaningless if the the parameter values deduced in this way fall 
outside the range of model’s validity. In some cases, these conditions 
of validity simply reflect restrictions on parameter values justifying 
certain mathematical approximations made to express a certain result 
in closed analytical form. Checking whether the various parameters 
deduced from data fitting do indeed satisfy the corresponding (math¬ 
ematical) restrictions or are consistent with ab initio estimated values 
is absolutely necessary. Nevertheless, not rarely current data analysis 
misses this important consistency step, which may be related to an 
insufficient discussion in the literature of the physical background of 
the various phenomenological models and the limited validity of the 
pertaining analytical formulas. Exposing the limits of applicability 
of such models frequently used in molecular electronics represents 
another aim of the present paper. In doing that, this paper is intended 
as an effort to make the community more aware of these limitations. 
It should by no means be understood as challenging the utility of 
model-based studies in gaining valuable conceptual physical insight 
into charge transport. 

Throughout, we will consider the case of zero temperature, as ap¬ 
propriate for off-resonant tunneling and symmetric coupling to elec¬ 
trodes (equal width parameters F/, = Tr = T due to “left” and “right” 
electrodes). In agreement with the fact that conductances of typi¬ 
cal molecular junctions are much smaller than the quantum conduc¬ 
tance G <C Go = 2e^//i = 77.48)iS, we will assume throughout T’s 
much smaller than relevant (molecular orbital) energy offsets relative 
electrodes’ Eermi energy eg = — Ef (Ef = 0). We will further 

assume eg > 0, because, for tunneling barriers, the corresponding 
formulas are simpler to write for electrons (n-type conduction); in all 
formulas presented below, eg should be understood as |eg| in cases 


of holes (p-type conduction) tunneling across negative barriers, as 
easily obtained by performing the charge conjugation transformation 
(eg = Eg — Eg —> Eg — Eg = —eg, e —> —e, m —>■ —m, see, e.g., ref. 

Likewise, since the tight-binding Hamiltonians of eqn (SH and 
(S3) are invariant under particle-hole conjugation (e.g. ref. [l9| and 
citations therein), the results for +eg and — eg coincide. 

While discussing in general the aspects delineated above, we will 
examine the benchmark case of conducting probe atomic force mi¬ 
croscope (CP-AFM) molecular junctions based on oligophenylene 
dithiols (OPDsjil as a specific example. 


2 Physical properties envisaged 


In this paper, we will mainly consider two physical properties (zero- 
bias conductance G and transition voltage Vj— ) and focus our atten¬ 
tion on how these properties vary with the molecular length d or size 
N across homologous series of molecules described within schematic 
models. Typical for nonresonant tunneling transport is a dependence 

G oc exp(—jSA) = exp (2) 


where /3 (j8) characterizes the exponential decay of the zero-bias con¬ 
ductance G with the molecular size (length). 

“Historically”,— the transition voltage V) was introduced as the 
bias at the minimum of the Fowler-Nordheim quantity log(//V^); 
the initial claim was that of a mechanistic transition from direct 
tunneling across a trapezoidal energy barrier {V < V)) to Eowler- 
Nordheim (field-emission) tunneling across a triangular barrier (V > 
V)).— As already noted in previous studies,i2i22i2i such a Fowler- 
Nordheim transition does not occur in molecular junctions. Physi¬ 
cally, Vt has nothing to do with the Fowler-Nordheim tunneling the¬ 
ory In particular, in molecular junctions the quantity log (//T^) 

does not linearly decrease with \/V at higher biases; this was a re¬ 
sult deduced for extraction of electrons from cold metals by intense 
electric fields— 

Still, transport measurements in molecular junctions do yield 
curves for log(//V^) exhibiting a minimum. Mathematically, the bias 
V = Vt at the minimum of log (//T^) coincides with the bias where 
the differential conductance is two times large than the nominal con¬ 
ductance—^ 


dl 

dV 


v=v, 


I 

= ^V 


v=v, 


(3) 


Eqn l|3 can be taken as alternative mathematical definition of Vt re¬ 
vealing at the same time its physical meaning. Using the minimum 
of log(//V^) plotted a function of V (or 1/V)il merely represents 
a mathematical trick to extract the V)-value of eqn (O- To eliminate 
the confusion that continues to exist in the literature on a mechanistic 
(Fowler-Nordheim) transition occurring at V = Vi, instead of Fowler- 
Nordheim diagrams log(//V'^) vs. V, it might be more appropriate 
to use diagrams V^/I vs. V, which have maxima exactly at the same 
V = Vt given by eqn so, instead of “transition voltage spectra” 

(“TVS”, Fig.[T]5), one can speak of a “peak voltage spectra” (“PVS”, 
Fig.IIJ;). 

As expressed by eqn (0, Vt represents a genuine nonequilib- 
rium property characterizing the charge transport through a molec¬ 
ular junction out of equilibrium. So, V) is qualitatively different from 
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Fig. 1 (a) Fitting the red I — V curve measured for a CP-AFM Ag/OPDl/Ag molecular junctioup^ with a cubic polynomial, eqn (|Tll, depicted 
by the blue line may appear to be quite satisfactory, (b) Nevertheless, for the same curves, the difference in the positions of the minimum of 
the Fowler-Nordheim (FN) quantity log (//F^) defining the transition voltage F, is substantially larger than experimental uncertainties 
{8Vt <0.1 V^). In the case illustrated here, the transition voltage obtained via eqn Q is even beyond the bias range accessed experimentally, 
(c) As alternative to the FN plot, one can inspect a “peak voltage spectrum”, namely the quantity F^// plotted vs. which exhibits a 
maximum located at F = Vf. (We show only the range F > 0 because the measured / — F curve is symmetric to a very good approximation.) 


the zero-bias conductance G; in fact, the latter (as well as other prop¬ 
erties commonly targeted in measurements, e.g., thermopower See- 
beck coefficient) can be expressed by properties of a device (let it be 
molecular junction or else) at equilibrium via fluctuation-dissipation 
theorem.— Interestingly, as revealed by recent studies^^ and sup¬ 
ported by measurements comprising thousands of molecular junc- 
tion s it is even more justified to consider F (rather than G) as a 
junction’s characteristic property; in a given class of molecular junc¬ 
tions, F( is much less affected by inherent stochastic fluctuations than 
individual I — V traces or low bias conductances G.— The aforemen¬ 
tioned should be taken as motivation why F/ represents a quantity on 
which the present transport study focuses. 
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Notice that the expansion coefficients can be expressed as 


(5) 

( 6 ) 


3 Models and relevant analytical formulas for 
symmetric I — V curves 

In this section we will examine a series of models widely employed 
in studies on charge transport through molecular devices. For some 
of these models, the dependence I = /(F) or other useful formulas 
can be given in closed analytical form. Worthy to be remembered, 
describing idealized situations, these analytical results hold only if 
certain conditions are satisfied. Whenever the case, the conditions 
of validity will be also given below. Other models examined in this 
paper rely on tight-binding Hamiltonians, for which the exact depen¬ 
dence / = /(F) can be obtained without substantial numerical effort. 

For reasons presented later, the expansion coefficients C 2 and C 4 
of the current up to the fifth order, eqn will also be provided for 
each case 

l{V) = Gv[\+C2V^ + c^V^'^+ff{y''^ (4) 

They can be directly inserted in eqn 0 and to compute the tran¬ 
sition voltages F ,3 and F (5 within the third and fifth order expansion. 


C 2 = . C 4 = I/F44 


where F 2 and F 4 have dimensions of voltages and allows one to ex¬ 
press the expansion in terms of dimensionless voltages 


/(F) = GF [1 -f (F/F2)^ -f (F/F4)'^] -F ff (f^) 


Next we briefly describe the models considered here and, when¬ 
ever possible, give the relevant analytical formulas for the low bias 
conductance G (Gq = 2e^ jh = 77.48/iS is the conductance quan¬ 
tum), the expansion coefficients C 2 and C 4 entering eqn (0 as well as 
the transition voltage Vt along with the conditions for their validity. 

(i) For biases e|F| < e^, within the Simmons WKB-based approx¬ 
imation for electron tunneling across a tunneling barrier of effective 
height Eg without lateral constriction (Simmons’s model), the current 
/ is given by eqn 0.'4 The Taylor expansion of the RHS of eqn 0 
allows one to deduce formulas for the zero-bias conductance G and 
the coefficients C 2 and C 4 entering eqn 0 . They are expressed by 
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eqn i), Gl, and {nil, respectively 


/(V) = - 
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- ( EiJ + ^ ) exp 11 -z^/ 1 + — 


eV 
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-pd). 


2eb 

G sa constant e^^'' 


= 0.0104167 


z (z^ — 3z — 3) 


C 4 = 0.0000325521 
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Validity condition: z = j3rf~4 


Above, £/ stands for the junction’s transverse area and 


a = 2 


\/2m* 


I ttl ^ 9 —] /2 I ffl 

= ao\ — = 1.025eV^'A — 


m 


mo 


(7) 

( 8 ) 

(9) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 


Oq = 2^y2mo/h, mo and m* being the Sommerfeld’s constant—, free 
and effective electron mass, respectively. Eqn dlU (like eqn l ll 8 t be¬ 
low) defines a dimensionless junction width z. The transition voltage 
for the Simmons and Simmons-based models has been investigated 
in a series of recent works. 

Eqn ll7]l and ([ 8 l( represent mathematical approximate results; in ad¬ 
dition to the restriction on the bias range specified, they only hold for 
barriers satisfying eqn d, requiring tunneling barriers sufficiently 
high and wide.— In fact, the applicability of Simmons’ approxima¬ 
tion is more restrictive than required by eqn d, as revealed by a 
more detailed analysis, 

(ii) Eor electron tunneling across a tunneling barrier with lateral 
constriction at biases e\V\li£B, the current I can be expressed by 
eqn d,— which applies provided that eqn ( 118b is satisfied. In ad¬ 
dition to the quantities G, C 2 and C 4 obtained by expanding the RHS 
of eqn d, the transition voltage V, can also be expressed in closed 
analytical form; see eqn d, d, and d 
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(14) 

(15) 

(16) 
(17) 


Validity condition: z = /3d~8 


(18) 


Noteworthy, the restriction imposed by eqn d to the results of 
eqn d and ( fTSl l is more severe than that expressed by eqn ( fT2t . 
which refers to a situation less appropriate for molecular electronic 
devices. Like eqn d, eqn d expresses the fact that the above 
results for model (ii) hold for physical situations where tunneling 
barriers are sufficiently high and wide. 


(iii) The highly off-resonant sequential tunneling— (superex¬ 
change limit—) across a molecular bridge consisting of a wire with 
N sites and one energy level (e^) per site represents a limiting case of 
the situation depicted in Fig.|2^. Provided that eqn d22b is satisfied, 
the current is expressed by eqn d; eqn d and d follow via 
straightforward series expansion. 
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( 21 ) 


Validity condition: EbI^ tf,,T 


( 22 ) 


Here, f/, is the hopping integral between nearest neighboring sites. 
Eqn ( 119b and ( 120b emerge from the Landauer formula, eqn d25b . using 
the transmission given in the RHS of eqn (S2) for arbitrary N. Notice 
that the power dependent transmission expressed by the latter yields 
a strict proportionality— 


eVt = 2ut£B £b (23) 

where the dimensionless quantity ut is the solution of the parameter- 
free algebraic eqn d, which results by inserting eqn d into 
eqn 13 

[N — 1/2) Ut 1^(1 + “f)^^ + (1 ~ 

= (1 — M()(l-fnf)^^ — (1 + M/)(1 — Mf)^^ (24) 

Eqn d22b expresses the physical fact that the above results for model 
(iii) hold for situations corresponding to strongly off-resonant tunnel¬ 
ing. 

(iv) A molecular chain consisting of N monomers, each monomer 
being characterized by a single orbital is schematically depicted 
in Fig. [^. Model (iii) presented above represents a limiting case 
(cf. eqn d22b ) of this model. To illustrate this, we give in eqn (S2) of 
the ESIf the transmission function computed exactly and within the 
sequential tunneling approximation for the case N = 2 Adapted to 
oligophenylene chains, £b should be taken as a model for the HOMO 
energy of an isolated phenylene unit, and tk as the coupling between 
HOMO levels of adjacent rings. (This is the motivation for choosing 
the subscript h.) 

Under applied bias V, two limiting cases can be considered. One 
limit is that of an applied field completely screened out by delocal¬ 
ized electrons (“metallic” molecule, Fig.|2};); the potential is constant 
across the molecule (V/ = 0), and the site energies 4 the same 
4 = £b- Another limit (no screening, “insulating” molecule) is that 
of a potential Vj- varying linearly from site to site between the values 
-|-V/2 and — V/2 at the two electrodes (Fig.[^), yielding r-dependent 
on-site energies — eV/. 
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(v) Molecular wires consisting of six-site rings with a single (7t- 
electron) energy level per site (CH-unit) depicted in Fig.|^ represent 
a tight-binding description of oligophenylene molecules. The cor¬ 
responding second-quantized Hamiltonian is given in the ESIf. As 
illustrated by Fig.l^t, li and t stand for intra- and inter-ring hopping 
integrals, respectively. 


a 



r=l r=2 r=3 r=4 


r=N 



i i 



t. 


C 



Fig. 2 Panels (a) and (b) represent schematically the tight binding 
models described by eqn (SI) and (S3). Panels (c) and (d) depict the 
cases of a constant potential and a linear potential drop, respectively, 
as limiting cases of spatial potential profiles between electrodes 
under a bias V. 


Although for models (iv) and (v) the dependence / = I(V) can¬ 
not be expressed in closed analytical form, it can easily be obtained 
numerically within the Landauer framework^ 

Go 

I[V) = — / T{E)dE (25) 

e i-eV /2 


Analytical formulas for the transmission function T{E) are lengthy 
even for F much smaller than Eb and f/,’s,^ 2 dl so they are of little 
interest in the present context. The numerical effort to compute T(E) 
via the Landauer trace formula is insignificant: it consists of a numer¬ 
ical integration and a matrix inversion, which is needed to compute 
the retarded Green’s function.— 

(vi) We will also examine the case of a generalized exponential 
transmission 


rg(F) =exp 


lE-eef 

AS 


(26) 


Results for this case, including an approximate analytical formula for 
the current valid for biases not considerably exceeding V) and 2 > A 
that can be deduced resorting to a Stratton-like approximatio n 
for transmission are presented below. The conductance G of eqn l l29b 
and the coefficients C 2,4 of eqn (ED and ED can be straightforwardly 
obtained from the fifth-order expansion of the RHS of eqn ED 
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(32) 

(33) 


The particular case 5 = 2, which corresponds to a Gaussian trans¬ 
mission, has occasionally been studied in the literatur e ^4^ ; results 
for this case are presented in the ESIf. 

The relation given below holds in general (not only for the Gaus¬ 
sian transmission 5 = 2) 


yapprox ^ 0.797483FJ'’''” = 1.02006<'^'’''“ (34) 

Here, and yl^PP’°^ represent estimates obtained from eqn (Ell 

and (|6j using the approximate coefficients and given 

above instead of C 2 and C 4 , respectively. An expression similar to 
yapprox (jgen given previously.— Notice that above the superscript 
approx refers to physical situations of strongly off-resonant tunnel¬ 
ing, mathematically expressed by the inequality A<eB. 


4 Results for symmetric I —V curves using 
generic model parameter values 

Detailed results for current (/) and transition voltages computed ex¬ 
actly and within the third- and fifth-order expansions (Vt,exact = ^t, 
V), 3 , and F;_ 5 , respectively) are collected in Fig. El STillEl S2, S3, 
S4,[6l|7l S5 and [ 8 ] Below we will emphasize some main aspects 
of the results obtained for the corresponding numerical simulations 
done with generic parameter values. 
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Fig. 3 The transition voltage Vt^exact computed by using the “exact” current for the Simmons model [termed model (i) in the main text], 
eqn 0 (panel (a)), along with the approximate estimates V, 3 and 14.5 obtained by inserting the expansion coefficients C2 and C4 of eqn (Uni 
and GD into eqn 0 and 0. Panel (b) show that, while V, 3 significantly deviates from Vt^exact (typical experimental uncertainties in Vt do not 
exceed ~ 10%^i^), 14 5 agrees with Vt^exact within a few percents {z = jSrf is a dimensionless barrier width). 




Fig. 4 (a) The transition voltage Vt.exact computed from the current of eqn ([H corresponding to the superexchange mechanism across a 
molecular wire modeled as chain of sites having a single orbital of energy Eg [termed model (iii) above], along with the approximate estimates 
V],3 and 14 5 obtained by inserting the expansion coefficients C2 and C4 of eqn ED into eqn 0 and 0 (panel (a)). Panel (b) shows that 14^3 
deviates from Vt^exact by up to 73 % (much larger than typical experimental uncertainties in I4 of ~ 10 %^i^), while V) 5 agrees with Vt^exact 
within at most 14 %. Notice that, within the sequential tunneling approximation, Vr^exact, 14.3, l^/,5 do not depend on the inter-site hopping 

integral 1 /,. 


4.1 The need to go beyond the parabolic conductance ap¬ 
proximation 


As visible in the aforementioned figures, the estimate I4 5, based on 
the fifth-order expansion of eqn 0, represents a good approxima¬ 
tion for the transition voltage V/^exact = 14 computed exactly for each 
model considered. For reasonably broad ranges of the model param¬ 
eter values, the relative deviations |l4 5 —141 /V) usually fall within 
typical experimental errors (.$ 10 %).^i^ This does not apply to the 
(over)estimate I4 3, whose deviations from the exact values are con¬ 
siderably larger that experimental errors. As visible in Fig.| 4 ]and S8, 
14.3 could be almost two times larger that Vt^exact- This demonstrates 
that the cubic approximation for current (or parabolic approximation 
for conductance), eqn 0, can only be used for qualitative purposes. 


4.2 The spatial potential profile as possible source of neg¬ 
ative differential resistance 

The impact of the spatial potential profile across a molecular junction 
is another issue worth to mention. Presumably, realistic potential pro¬ 
files fall between the two situations depicted in panels (c) and (d) of 
Fig. [2 (a constant potential and a linearly varying potential, respec¬ 
tively); so examining the differences in the results obtained in these 
two cases may be taken to assess how important is to exactly know 
the actual potential profile in a real case. The results presented in 
Fig. S 2 , S 3 and S 4 indicate a weak effect at biases not much larger 
than Vt- Obviously, this is a further plus of transition voltage spec¬ 
troscopy. 

The manner in which the potential varies across a molecular junc¬ 
tion is important only at biases substantially higher than V). This is 
shown in Fig. S 2 , S 3 , S 4 and especially in Fig.[^ which, to the best 
of our knowledge, demonstrates a qualitatively new effect. As vis¬ 
ible there, instead of a current plateau occurring for a potential flat 
across the molecule beyond bias values at which molecular orbital 
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Fig. 5 Transition voltage Vt computed exactly for a molecular wire whose monomers are modeled as single sites characterized by a single 
orbital of energy eg- For arbitrary nearest-neighbor hopping integral values f/,, this is termed model (iv) in the main text, which becomes 
model (iii) in the limit f/, <C £b- Notice the almost insignificant change 14 —14 + 8Vi when the spatial potential profile across the chain 
changes from constant (Fig.[^) to a linear drop (Fig.[2jl). The differences between the dashed and solid lines in panel (a) represent the 
absolute changes 5l4. The relative changes 514/14 are shown in panel (b). 


energies (£ 1 . 2 ,...) become resonant to the Fermi level of one electrode 
(eVres = 2 £i. 2 ,„,)^, the current decreases for a linearly varying po¬ 
tential. We are not aware of a previous indication on such a nega¬ 
tive differential resistance (NDR) effect related to a specific potential 
variation across a molecular junction. Noteworthy, this is an NDR 
effect for uncorrelated transport computed within the limit of wide¬ 
band electrodes. Previous work on uncorrelated transport indicated 
finite electrode bandwidths and energy-dependent electrode density 
of states as possible sources of NDR.— 

4.3 Gaussian transmission versus Lorentzian transmis¬ 
sion 

The comparison between the cases of Lorentzian (obtained as a par¬ 
ticular case of model (iv) for N = I, discussed in detail previously, 
e.g., ref. 113 ) and Gaussian transmission (a particular case of eqn d26b ) 
reveals that, even systems wherein the transport is dominated by a 
single energy level, may have qualitatively different physical proper¬ 
ties. 

As visible in Fig.|7]and S5 and eqn (S7), V) strongly depends on the 
width parameter A and is inversely proportional to By contrast, 
in the case of Lorentzian transmission, V) linearly increases with Eg 
and is nearly independent of the width parameter F as long as F <C £b 
[ eV4 = 2Eg/V3 + ^(F/£B)^].i2i2ii22 In fact, both the proportional¬ 
ity to Eg and the insensitivity to F-variations for F ^ £g is a com¬ 
mon feature shared by the above models (iii-v), while a Vt nearly 
1 /2 

proportional to Eg at large Eg for models (i) and (ii) is the con¬ 
sequence of the fact that the corresponding transmissions are simi¬ 
lar to that of eqn ( 126b with 5 = 1/2 (cf. eqn ( 117b and ( 130b ). Com¬ 
mon features for these two types of transmissions are the inflection 
point of the I —V curves (“current steps”. Fig. [^) located at biases 
Vres = 2Eg /e where the level becomes resonant to electrodes’ Fermi 
energy, and the fact that the associated maximum in the differen¬ 
tial conductance has a width proportional to the transmission widths; 
SVres A (Fig.[^), similar to SVres F for Lorentzian transmission. 
The opposite dependence of V) on Eg for Lorentzian (V4 Eg) and 
Gaussian (Vt 1 /Eg) transmission is relevant for the correctly inter¬ 


preting gating effectsii^ or the impact of electrodes’ work func- 
tion .22 

Concerning the approximate estimates for Vt^exact, Fig.ITJt and S5b 
indicate a behavior similar to that found in the above cases: 1 / 5 rep¬ 
resents an accurate estimate for the exact Vt^exaci, while F ,,3 is quan¬ 
titatively unsatisfactory. 

5 Asymmetric I —V curves described within 
the Newns-Anderson model 

Let us briefly examine the case where, unlike those assumed through¬ 
out above, the I — V curves are not symmetric under bias reversal, i.e., 
I(—V) ^ ~I(V). It is clear that the asymmetry with respect to bias 
polarity reversal cannot be satisfied within the framework provided 
by eqn (0 and even powers in V should also be added 

/(F) = /3(F)-|-^(f‘*) = 

Gv(^\+ciV + C2V^'^+ff(v'^^ (35) 

/(F) = Is{V) + ff(v^'^=Gv(^l+ciV + C2V^ 

+ c2V^+C4V‘^'j+ff(v^'^ (36) 

By including only the lowest even power (i.e. ci 7 ^ 0, C 3 = 0)^, it is 
possible to account for an asymmetry /(F) —/(—F). However, do¬ 

ing this, i.e. using eqn HD. has an important drawback. By inserting 
eqn (ED in eqn m one easily gets 

V),3+ = -Vr, 3 - = (37) 

exfci 

i.e., transition voltages of equal magnitudes for positive and nega¬ 
tive bias polarities, and this holds no matter how large is ci. Note¬ 
worthy, the asymmetry /(F) 7 ^ —/(—F) of an / — F curve does not 
automatically imply Vt- 7 ^ —Vt-i-. Experimenta l and theoreti¬ 

cal— results show that situations wherein Vt- 7 ^ — Vt+ are physically 
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Fig. 6 Tunneling current computed for molecular wires [model (iv), eqn (SI)] for a constant (Fig. [2};) and linearly varying (Fig.|^) potential. 
The linear potential drop has an impact at higher biases more pronounced than that visible in Fig. S2, S3 and S4; displacements of the current 
step positions (which are located at on-resonance situations, eV = 2ey, in the case of a constant potential) and negative differential effects, e/s 
(j > 1) represent the absolute values of the highest occupied orbital energies measured from the electrodes’ Fermi energy. For a two-site 
chain, these values are £12 = (cf- inset of panel (a)). Notice the logarithmic scale on the y-axis. 




Fig. 7 Results obtained using the Gaussian transmission of eqn (S4) showing the dependence on Sb of the transition voltage computed exactly 
(Vt), and using its approximate estimates: Vt^appmx of eqn (S7) and and V, 5 obtained from eqn (|5}, (S 8 ) and (S9). Notice that 

both y/PP''°^ and V, 5 represent good approximations for large values of the ratio Eg/A. 


relevant; this clearly demonstrates that the “generic parabolic” de¬ 
pendence of the conductance of eqn d35b is an unsatisfactory approx¬ 
imation. 

To find the counterpart of eqn (HJl for the asymmetric case one 
has to solve an algebraic fourth order equation obtained by inserting 
eqn l l36b into eqn ([3]l 

34v,^ -f 2 c 3-f C 2 = 1 (38) 

This yields asymmetric transition voltages V), 5 -|- / Their an¬ 

alytical expressions are too long and will be omitted here. 

Starting from the models discussed above, asymmetric I — V 
curves and Vt,+ / —Vt - can be obtained by allowing a bias-induced 
energy shift 

eB-*eB{V) = eB + YeV (39) 

where —l/ 2 <y< 1/2 is a voltage division factor (see, e.g., ref. fl^ 
and citations therein). 


Although calculations are straightforward, being too lengthy, the 
counterpart of the formulas given for symmetric 1 — V curves in 
the preceding section will not be given here. Instead, we will re¬ 
strict ourselves to the Newns-Anderson (NA) model— ^ within the 
wide-band limit, which assumes a single level of energy £b(F) and 
Lorentzian transmission. For (off-resonant) situations of practical in¬ 
terest (F <C Eb, biases up to ~ l.SEg/e) 29,31,46,53- ^^ exact formulae 
for the current and Vi± have been deducedi^; see eqn (SIO) and (Sll) 
in the ESIf. The fifth-order expansion in eqn (S10) straightforwardly 
yields the expansion coefficients entering eqn I36b 

Cl = -2y ; C 2 = ^ -b 3y^ (40) 

C3 = -y(l+4y2j ; c4 = :j^+5y2Q-l-'rj (41) 

where Cj = cj (eb/c)^ (j = 1 to 4). By inserting the above expression 
in eqn ( I37l l and OSb . the approximate values Ff^ 3 ± and V) 53 - for both 
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Fig. 8 Results obtained for a Gaussian transmission showing that, while the voltages where the I — V curves (panel (a)) exhibit an inflection 
point (or maximum of the differential conductance dljdV of width proportional to A, panel (h)) have a dependence similar to the case of a 
Lorentzian transmission (namely e|y| = 2£g corresponding to the alignment of the energy level Eg with the Fermi level of the biased 
electrodes), the transition voltage decreases with increasing Eg, in contrast to the case of Lorentzian transmission, where it increases with Eg 
(eV, = 1.154eB)i^. 



Fig. 9 Results for a single level and Lorentzian transmission. Whether computed exactly, using the third- or fifth-order expansions eqn (SIO), 
( 135b . and l l36b . respectively) I — V curves are asymmetric with respect to the origin for a nonvanishing voltage division factor y (panel (a)). 
The transition voltages (minima of the Fowler-Nordheim quantity in panel (b)) computed exactly from eqn (Sll) for opposite polarities are of 
different magnitudes (V/+ 7 ^ The fifth-order expansion correctly describes this inequality (V(_ 5 + 7 ^ |V), 5 -1), while the third-order 

expansion incorrectly predicts V/ 3 + = IV'r.s-l. 


bias polarities can be obtained. The comparison with the exact 'Vt± 
obtained from eqn (S11) is depicted in Fig. S 6 andllOl 

To end this subsection, lower order expansions (eqn ( 135b and ( 136b ) 
for the asymmetric I — V curves considered above appear to be more 
problematic than for symmetric cases. Fig. shows that even the 
fifth-order expansion provides a satisfactory description for both bias 
polarities only for weak asymmetries (|y|~ 0 . 1 ). 

6 Case of molecular junctions based on 
oligophenylene dithiols 

Analytical expressions for current and conductance like those given 
above are very often utilized by experimentalists to fit their mea¬ 
surements. However, and this is one important point on which the 
present work wants to emphasize, these formulas are only valid un¬ 
der well defined conditions/restrictions, and utilization beyond these 
conditions makes no sense. To exemplify, in this section we will an¬ 
alyze recent transport data obtained for molecular junctions based on 


oligophenylenes using silver electrodes Ag/OPDs/Ag containing up 
to four phenyl rings (1 < A < 4).— 

Experimental studies of OPD-based molecular junctions with sil¬ 
ver electrodes containing 1 < A < 4 phenyl rings yielded a value of 
the conductance attenuation factor j3 ~ 1.56 per ring;— for a phenyl 
ring size d\ = 4.3 A, this corresponds to j3 ~ 0.36 A —. This j3- 
estimate can be invoked to immediately rule out the tunneling bar¬ 
rier picture underlying models (i) and (ii) discussed above as valid 
framework to analyze the charge transport in these junctions. Out 
of the phenylene dithiol species studied in ref. it is only the four- 
member species, d ^ d^ = 4rfi, that satisfies condition of eqn GU, 
pd = 6.24 > 4. In fact, as discussed in detail recentl y the Sim¬ 
mons model in its original formulation— does not account for the fact 
that, in molecular junctions, electron motion is laterally confined. 
When lateral constriction is incorporated into theory, the condition 
becomes more restrictive. Instead of eqn G3, one should apply 
eqn GD, which is invalidated even for A = 4. 

The presently employed parameter Eg represents an effective 
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Fig. 10 Results for a single level and Lorentzian transmission. Transition voltages for positive and negative biases computed exactly Vt± from 
eqn (SI 1), and using the the third- or fifth-order expansions [Vf.srp, eqn ( I35l l. and eqn 06l l. respectively]. Notice the incorrect prediction 
Vt 3 + = |Vf,3-1 - The fifth-order expansion provides good estimates — Vt± only for small 7 ; at larger 7 , only one polarity is satisfactorily 
described (e.g., Vf 5 _ Vt- for 7 > 0 ). 


barrier height, which also embodies image effects. According to 
eqn GD and to “explain” a certain /3-value adjusted to fit exper¬ 
imental data, both Eb and the the effective mass m* can empirically 
be “adjusted” to make theory to “agree” with experiment. However, 
what cannot be “manipulated” in this way is the //-value itself. The 
restriction expressed by eqn l ll 8 t represents the mathematical condi¬ 
tion for a valid description within the tunneling barrier model, and 

it cannot be modified whatever “appropriate” is the choice of Eg and 
* 

m . 

The argument based on the (too small) /3-value also demonstrates 
that a description based on the superexchange limit of tunneling, 
which underlines model (iii),— is impossible. For /3 = 1.56, eqn ll20t 
yields r/,/eg = 0.458, which is not much smaller than unity, as re¬ 
quired by eqn i22\ . Again, whatever “appropriate” the adjustment 
of the parameters Eb and f;,, they cannot be chosen to satisfy the 
condition required by theory for a valid superexchange mechanism, 
because the experimental data “fix” the value of /3 and thence via 
eqn iQj the ratio t^/EB- So, one should go beyond the superexchange 
limit. 

The most straightforward generalization is to use model (iv) for 
a chain with having one level per monomer, for which theory does 
not impose restrictions to the ratio between Eb and tf,. Within a mi¬ 
croscopic description based on model (iv), f;, represents the hopping 
integral coupling the two HOMOs of two adjacent benzene rings. It 
can be deduced from the difference between the ionization energies 
/i and I 2 of benzene and biphenyl: 2f;, = h— h- Quantum chem¬ 
ical calculations based on refined ab initio methods, EOM-CCSD 
(equation-of-motion coupled clusters singles and doubles), as im¬ 
plemented in CFOUR— and OVGF (outer valence Green’s func¬ 
tions)^^ done with GAUSSIAN 09^ allowed us to reliably de¬ 
termine /i and 12- The values are /i = 9.219eV and 9.197 eV, and 
I 2 = 8.217 eV and 8.204eV, respectively. They allow to compute 
the hopping integral needed for calculations based on model (iv): 
tij = 0.501 eV and 0.497 eV, respectively. The agreement between 
these two t/,-values demonstrates the reliability of the ab initio es¬ 
timates. The direct ab initio determination of Eb (relative align¬ 
ment between HOMO and electrodes’ Fermi level) is an issue too 
difficult to be addressed here. It can be determined from the ex¬ 
perimental value ci 1.15 V for Ag/OPDl/Ag.^i21 The equa¬ 


tion Eb = eVtV2/2 — yields Eb = 0.996 eV. Exact currents for the 
Hamiltonian of eqn (SI) can be computed numerically, and this al¬ 
lows to obtain V) for various N’s. Results based on model (iv) us¬ 
ing this ab initio estimated 1/, = 0.5 eV and adjusting Eb to the value 
Eb = 0.996 eV fixed by the experimental value — 1.15 V— are 

completely unacceptable; not even an exponential decay of the con¬ 
ductance with increasing size N can be obtained. 

In view of this disagreement, we have attempted to keep only one 
of the above parameters (either q, = 0.5 eV or Eb = 0.996 eV) and 
to determine the other by fitting the experimental value /3 = 1.56.— 
Fixing Eb = 0.996 eV yields tfj = 0.387 eV, while fixing ?/, = 0.5 eV 
requires a value Eb = 1.288 eV. The results obtained in this way are 
depicted in Fig. S 6 andll II respectively, and they show that none of 
these empirical changes can make model (iv) to provide a satisfac¬ 
tory description. Differences between the cases of a flat potential and 
a linearly varying potential are insignificant (Fig. S 6 c and fTTll: so, 
the lack of information of the real spatial potential profile across junc¬ 
tions cannot be advocated as possible source of this disagreement. 

A description based on model (v), which considers the phenyl 
rings explicitly, represents the highest reasonable refinement of a 
tight-binding approach to charge transport in OPDs junctions. To 
determine, f/, which within this model is a characteristic of a ring, 
we use the lowest singlet excitation energy Egxc\ it can be expressed 
as Egxc = 2r;. Using aug-cc-pVDZ basis sets, we have determined 
Eexc from EOM-CCSI>S and SAC-CI (symmetry adapted clus¬ 
ter/configuration interactionjiS calculations; the values are E^xc = 
5.154 eV and Egxc = 5.051 eV, respectively. They are in good agree¬ 
ment with the experimental values: 4.9e V ^^ 1 ^^ and 5.0 eV— So, we 
use the value q = 2.5 eV, which also agrees with the overall descrip¬ 
tion of the excitation spectrum of benzene as well as with polyacety- 
lene data—. To determine t, we use again the difference between 
the lowest ionization energies of benzene and biphenyl. (Notice that 
the parameter t of model (v), which is the hopping integral between 
neighboring C sites belonging to adjacent phenyl rings (Fig.| 2 j 3 ), is 
different from the parameter f;, of model (iv), which is the hopping 
integral between HOMOs of adjacent phenyl rings (Fig.[^).) To re¬ 
produce I 1 —I 2 — 1 eV (see above), a value t = 3.677 eV is needed. 
The very fact that t = 'i.&ll eV is larger than q = 2.5 eV is unphysi¬ 
cal; the hopping integral t associated with a single C-C bond cannot 
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Fig. 11 Results for model (iv) using the ab initio value 4 = 0.5 eV and the value eg = 1.288 eV obtained by fitting the experimental 
conductance tunneling attenuation coefficient j8 = 1.56 (panel (a))^ The agreement between the calculated Revalues (panel (b)) and 
experiment is poor. The spatial potential profile across these junctions do not notably affect Vj, as illustrated by the results obtained for 
constant and linearly varying potentials (panel (c)). The points (panels (b) and (c)) and error bars (panel (b)) represent experimental results for 
CP-AFM Ag/OPDs/Ag junctions.^ 


be larger than the hopping integral f,- associated to neighboring carbon 
atoms in an aromatic ring characterized by C-C distances shorter than 
of a single C-C bond (Fig.[ 2 j 3 ). It is a clear expression that, param¬ 
eters of tight-binding models (however refined they are) cannot be 
adjusted to satisfactorily oligophenylene molecules. This finding is 
in line with recent studies demonstrating limitations of tight-binding 
approaches for molecules of interest for molecular electronics.— 

Using the ab initio values t = 3.677 eV and f/ = 2.5 eV and adjust¬ 
ing Eb to fit ~ 1.15 V yields eg = 3.433 eV. This parameter set 

is not even able to qualitatively describe the exponential decay with 
N of the conductance. Therefore, we have also considered two al¬ 
ternative parameter sets: eg = 6.74 eV, t; = 2.5 eV, and t = "i-Gll eV 
(Fig.nil and eg = 3.433 eV, q = 2.5 eV, and t = 0.964eV (Fig. S7). 
To get the first set, instead of fitting V) ~ 1.15 V (as done above), 
we have adjusted eg by imposing jS = 1.56. To obtain the second set, 
we have adjusted eg and t to fit the experimental values jS = 1.56 and 
~ 1.15 V. Notice that the value t = 0.964 eV is much smaller 
not only than the ab initio estimate t = 3.677 eV and also smaller than 
t, = 2.5 eV, which cannot be explained by the fact that a C-C single 
bond is slightly shorter than an aromatic one— As visible in Fig.ll2l 
and S7, none of these two sets provides a satisfactory description. 

The (in)accuracy of Vj-estimates based on the third- and fifth-order 
expansions of eqn 0 and @ using parameter values specific for 
OPDs junctions, which is illustrated in Fig. S 8 , is comparable to that 
encountered in the previous situations (Fig. [3lSl,El|7]and S5). 

To conclude this part, neither representing OPDs junctions as 
chains of phenyl rings merely considering coupled HOMO’s of 
neighboring rings within model (iv) nor a full treatment of the cou¬ 
pled rings at tight binding level is able to provide an acceptable de¬ 
scription of the transport experiments.^ Limitations of tight binding 
descriptions of molecules of interest for molecular electronics have 
been recently pointed out.— 

The data collected in [T] demonstrate why a description based on 
a single level and Gaussian transmission should also be ruled out 
for OPDs junctions. The difference between the values eg deduced 
from transport data— and from ultraviolet photoelectron spectroscopy 
(UPS>2i for Ag/OPDl/Ag is enormous in the only case {N = 1) 
where the latter is available. In addition, eg, which represents the 


N 


A 


expt 

1 

1.15 

1.273 

3.061 

1.1 b 

2 

1.00 

1.265 

3.378 


3 

0.84 

1.212 

3.615 


4 

0.70 

1.118 

3.639 



Table 1 Results obtained by modeling Ag/OPDs/Ag by assuming a 
Gaussian transmission, eqn (S4). experimental values^; result 
of ultraviolet photoelectron spectroscopy (UPS)^. 


HOMO energy offset relative to the Fermi level, is found to increase 
with increasing N, which is completely unphysical.— 

7 Additional remarks 

To avoid possible misunderstandings related to the present paper, a 
few remarks are in order, however. 

The above considerations refer to existing models, which are cur¬ 
rently used for interpret experimental data for the charge transport 
via tunneling in molecular junctions. Our main aim was to present 
a list of formulae that can be used for data processing along with 
the pertaining applicability conditions. This paper is intended as a 
working instrument enabling to check whether certain experimental 
transport data are compatible or not with one of the existing mod¬ 
els. By presenting the benchmark case of OPD-based junctions, we 
mainly aimed at illustrating how experimental data could/should be 
utilized to (in)validate a certain theoretical model, not merely check¬ 
ing whether I — V measured curves can be fitted without checking 
whether the values of the fitting parameters are acceptable. An ex¬ 
haustive analysis of transport data available in the literature for var¬ 
ious molecular junctions within the presently considered models is 
beyond the present aim; this may make the object of a (review) paper 
of interest on its own. 

The presently considered models disregard important physical ef¬ 
fects (e.g., details of interface molecule-electrode couplings or metal- 
induced gap states). It may be entirely possible that none of these 
models is entirely satisfactory just because such effects are signif- 
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Fig. 12 Results for the conductance (panel (a)) and transition voltage (panel (b)) computed within model (v) using the values r, = 2.5 eV and 
1 = 3.677 eV deduced ah initio and adjusting Eb = 6.74 to fit the tunneling attenuation factor /3 = 1.56. In panel (b), the points and error bars 
represent experimental results for CP-AFM Ag/OPDs/Ag junctions.— 


leant. A corresponding extension of the theoretical models to pro¬ 
vide experimentalists with simple formulae enabling data process¬ 
ing/interpretation, obviating demanding microscopic transport calcu¬ 
lations remains a desirable task for the future. 

In this paper we did not discuss the (off-resonant limit (F ^ Eb) 
of the) Newns-Anderson model in much detail. Relevant formu¬ 
lae (eqn (SIO) and (Sll)) are given in the ESIf. As shown in re¬ 
cent studie s cases where this description applies and is jus¬ 

tified microscopically exist. There, the transition voltage can be 
used to directly extract the orbital energy offset using a formula 
(eVt = 2eb / t/S = 1.155 eVf for y = 0i2) close to that initially claimed 
(eVt = Eb)^- However, the results presented for the various models 
analyzed above have demonstrated that V) cannot be uncritically used 
to straightforwardly deduce the energy alignment of the dominant 
molecular orbital. This does not diminish the importance of Vf. As 
emphasized in Section V) is an important property characterizing 
the nonlinear transport. By studying its behavior across homologous 
molecular classes (i.e., varying the molecular size d or n) or under 
gating (i.e., varying Eb), the (in)applicability of a certain model can 
be concluded. Moreover, Vt turned out to be a key quantity, as it 
allowed to reveal that charge transport across different experimental 
platforms and different molecular species exhibit a universal behav¬ 
ior, which can be even formulated as a law of corresponding states 
(LCS) free of any empirical parametersiii. 

In order to demonstrate that the cubic polynomial approximation 
is insufficiently accurate to describe V), we have presented in Fig. a 
a raw I — V trace measured on a CP-AFM junction^. Making the 
point in connection with Fig.[T]is only possible by using neat, smooth 
curves; typical raw I —V curves measured for single-molecule junc¬ 
tions are blurry and therefore inadequate for this purpose. A compar¬ 
ison between junctions consisting of a single molecule and a bundle 
of molecules is of interest on its own; still, we note that, as long as 
the charge transfer occurs via off-resonant tunneling, the similarity of 
the transport properties of these two types of junctions is expressed 
by the aforementioned LCS—. Considering transport through CP- 
AFM junctions within the models discussed above amounts to ne¬ 
glect proximity effects due to other (identical) molecules in the bun¬ 
dle. This assumption may certainly fail in some cases. However, 
the fact that OPD-based CP-AFM junctions obey this LCS^ can be 


taken as a strong indication that the failure of the various models in 
case of these molecular devices discussed in Section is not related 
to the approximate description in terms of of a bundle of independent 
molecules. 

8 Conclusion 

With the manifest aim of providing experimental colleagues a com¬ 
prehensive working framework enabling them to process and inter¬ 
pret measurements of transport by tunneling in molecular junctions, 
this paper have presented a detailed collection of analytical formu¬ 
lae, emphasizing on the fact (less discussed in the literature) that 
these formulae only hold if specific conditions of applicability are 
satisfied, which often impose severe restrictions on the model pa¬ 
rameters. From a more general perspective, the theoretical results 
reported above have demonstrated that: 

(i) The often accepted idea of a generic parabolic V-dependence 
of the conductance as fingerprint of transport via tunneling emerged 
from studies based on high and wide energy barriers. This picture 
misses a microscopic foundation for tunneling across molecules char¬ 
acterized by discrete energy levels. At biases of experimental inter¬ 
est, within all the models examined here, the description based on a 
third-order expansion 1 = /(V) of transport measurements in molec¬ 
ular junctions is insufficient for quantitative purposes. To illustrate, 
we have shown this by systematically analyzing the fifth-order expan¬ 
sions, which turned out to be reasonably accurate at least for biases 
up to the transition voltages. 

(ii) Merely adjusting model parameters (e.g., within a tunnel¬ 
ing barrier picture by claiming renormalization effects due to im¬ 
age charges or effective mass) does not suffice to “made” a model 
valid for describing a specific molecular electronic device. On one 
side, models are normally valid only under specific conditions, which 
these parameters should satisfy. On the other side, the parameter val¬ 
ues deduced from fitting experimental data should be consistent with 
ab initio estimates. Parameters for tight-binding models can be eas¬ 
ily estimated via reliable ab initio quantum chemical calculations, as 
illustrated in this study. 

(iii) The experimentally measured values of the /3-tunneling coef¬ 
ficient of the low bias resistance can be inferred for quickly assessing 
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the inapplicability of certain tunneling models. This turned out to 
be the case for OPDs junctions, whose (too small) value (jS = 1.56 
or jS = 0.36A ) deduced from experiments is incompatible with 

descriptions based on tunneling barrier or superexchange mecha¬ 
nism. Because molecular junctions based on other aromatic aromatic 
species have often even smaller /3’s (jS ~ 0.2 A ' descriptions 

based on those models should be excluded. In such cases, uncritical 
application of the mathematical formula of I vs. V is meaningless, 
even if the measured / — V curves can be satisfactorily htted. 

(iv) For biases not too much higher than the transition voltage, a 
realistic description of the spatial potential profile across a molecular 
junction appears considerably less important that usually claimed. 
This is no longer the case at high biases, where the spatial poten¬ 
tial prohle may be responsible for qualitatively new phenomena, e.g., 
negative differential resistance (cf. Fig. [ 6 ). 

(v) The fact that the estimate Vj 5 ~ Vt based on the fifth-order 

/(F)-expansion turned out to be acceptable in many cases can be 
of practical help; eqn 0 (or eqn J36l l) can be employed to process 
noisy experimental I—V curves, for which straightforwardly redraw¬ 
ing measurements as \og{IjV'^) vs. V (or V^// vs. V) may yield sub¬ 
stantial uncertainties to determine the minimum (or maximum) posi¬ 
tion. Fitting I — V curves with fifth-order polynomials can be used 
to extract the coefficients £ 2 , 3,4 and to estimate Vj ~ ^ by means of 

eqn 0 or l IMt . 

(vi) Finally, we refer to situations where a successful description 
based on a single level and Lorentzian transmission (also known as 
the Newns-Anderson model) has been concluded^i^2i2£. The fact that 
the I—V curves could be very well fitted within the Newns-Anderson 
framework was not the only argument leading to that conclusion; 
equally important was that the energy offset Eb deducing from fitting 
the I—V data has been correlated with ah initio estimates-^i^i and/or 
independent experimental information.— Even if the transport would 
have been dominated by a single level, the conclusion of those stud¬ 
ies would have not emerged in case of, e.g. a Gaussian transmission, 
because of the completely different dependence of Vt on eg. 

With the (few) examples mentioned in the preceding paragraph, 
we want to end by reiterating an idea already presented in Introduc¬ 
tion; while attempting to make the community more aware of the lim¬ 
its of applicability of the various models utilized, the present paper 
did by no means intend to challenge the overall usefulness of model- 
based studies in gaining conceptual insight into the charge transport 
by tunneling at nanoscale. 
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